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I.  INTRODUCTION 


Interferometers  measure  the  correlation  of  radiation  fields  that  have  taken 
different  paths  through  space.  Astronomers  have  widely  used  interferometers  because  the 
angular  resolution  of  the  device  goes  as  k/B.  where  X  is  the  wavelength  of  the  radiation 
and  B  is  the  separation,  or  baseline,  of  the  detector  elements  in  the  interferometer.  Thus, 
for  a  given  wavelength  it  is  possible  to  improve  an  interferometer’s  resolution  by 
increasing  its  baseline,  and  though  there  are  technical  difficulties  associated  with  this,  it 
is  nonetheless  an  effective  means  of  improving  resolution  as  opposed  to  increasing  the 
diameter  of  a  single- aperture  imaging  system. 

Interferometers  used  by  astronomers  come  in  two  broad  classes:  amplitude 
interferometers  and  intensity  interferometers.  Each  can,  within  its  own  set  of  limitations, 
be  used  to  image  remote  objects.  This  report,  which  was  funded  as  a  Central  Research 
Project  by  the  Science  and  Technology  Division  of  the  Institute  for  Defense  Analyses, 
gives  an  introduction  to  the  subjects  of  amplitude  and  intensity  interferometry. 

Obviously,  there  are  many  references  that  discuss  interferometry.  I  have  drawn 
from  a  few  to  compile  this  report.  Born  and  Wolfs  classic  text  on  optics  is  a  first-rate, 
broad  reference  to  the  field  [1].  Boden’s  chapter  in  the  1999  Michelson  Summer  School 
course  notes  is  an  excellent  introduction  to  amplitude  interferometry  [2].  The  original 
article  on  the  subject  of  intensity  interferometry  by  Hanbury  Brown  and  Twiss  provides 
unique  intuition  into  the  subject,  and  it  has  historical  significance  [3].  Thompson,  Moran, 
and  Swenson’s  text  on  radio  interferometry  in  astronomy  is  a  single,  comprehensive 
reference  on  the  various  interferometric  techniques  used  by  astronomers  [4]. 

The  response  of  an  amplitude  interferometer  and  its  application  as  an  imaging 
device  is  given  in  Section  II  of  this  report.  Section  III  presents  the  response  of  the 
intensity  interferometer  and  shows  how,  as  an  imaging  tool,  it  is  limited  to  measuring 
source  sizes  in  a  model-dependent  manner.  Extensions  to  intensity  interferometry  are  also 
introduced  in  the  chapter,  in  particular  the  methods  of  Gamo  [5]  and  Mehta  [6]  and  the 
experiment  by  Beard  [7],  all  of  which  superimpose  a  known  reference  beam  on  a  source, 
enabling  intensity  interferometers  to  fully  image  a  source  distribution. 
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The  interferometric  imaging  techniques  covered  in  this  document  have  the 
functionality  of  a  lens  and  focal  plane  array  when  imaging  distant  objects.  They  further 
have  the  flexibility  to  improve  resolution  by  expanding  the  baseline  of  the  interferometer 
beyond  the  size  scales  of  traditional  lens/focal  plane  array  systems.  Interferometers 
consequently  have  great  potential  when  high-resolution  images  of  distant  objects  are 
desired. 
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II.  AMPLITUDE  INTERFEROMETRY 


Consider  the  interferometer  pictured  in  Figure  1 .  The  radiation  field  from  a  source 
falls  on  two  apertures,  A;  and  A?,  separated  by  baseline  B.  The  fields  that  enter  Ay  and  A? 
are  combined  after  passing  through  the  two  arms  of  the  interferometer,  which  contain 
delays  dj  and  ch.  The  detected  intensity  is  the  square  of  the  sum  of  the  field  amplitudes  at 
the  point  that  they  are  combined.  The  relative  phase  of  the  fields  at  that  point  will 
determine  if  the  field  amplitudes  constructively  or  destructively  interfere. 

Such  an  interferometer  could  be  constructed  with  a  pair  of  linear  detectors* 
followed  by  cables  of  lengths  di  and  cf,  which  are  fed  into  circuitry  whose  power  output 
is  the  square  of  the  sum  of  the  signal  amplitudes  at  the  circuit’s  input.  The  device  could 
also  be  a  Young’s  double-slit  interferometer  where  the  two  apertures  are  cut  in  an 
aperture  screen  and  the  field  combination  is  seen  on  an  image  screen  positioned  behind 
the  apertures.  In  this  case,  the  relative  delays  di  and  cf  would  be  represented  by  the 
position  viewed  on  the  image  screen  along  the  line  parallel  to  the  apertures’  baseline. 
Whatever  its  manifestation,  this  class  of  device  is  referred  to  as  an  amplitude 
interferometer,  because  the  detected  intensity  results  from  the  square  of  the  sum  of  the 
field  amplitudes  detected  by  the  interferometer. 

A.  RESPONSE  TO  A  DISTANT  POINT  SOURCE 

Consider  a  quasi-monochromatic  point  source  with  wave  number  k  =  2n/X  (where 
X  is  the  wavelength  of  the  source)  at  an  angle  9  from  the  normal  to  the  baseline  of  the 
device.  Assume  that  9  is  small  (sin 9  «  9),  that  the  source  is  distant  so  that  its  wave-front 
is  a  plane  wave  at  the  interferometer,  and  that  the  apertures  and  responses  in  each 
interferometer  arm  are  equal.  The  amplitude  of  the  combined  signal,  anet ,  is  determined 
by  the  relative  phase  of  the  signals  in  the  spectrometer  arms,  a]  and  a 2  (common  phase 
factors  can  be  absorbed  into  an  overall  constant  and  are  of  no  consequence): 

a„et  =  1  +  a2  ~  cos(£Z?9  +  kdx  -  cot)  +  cos (kd2  -  cot)  . 


A  linear  detector  is  one  whose  voltage  output  is  proportional  to  the  amplitude  of  the  input  field. 


In  an  amplitude  interferometer,  we  will  observe  the  tzme-averaged  power,  /,  of  the 
combined  field,  where  the  time  average  is  over  a  period  T  that  is  much  larger  than  the 
period  of  the  radiation.  Note  that  this  is  not  the  instantaneous  power,  which,  for 
frequencies  above  about  1  GHz,  is  not  practical  to  observe.  Since  /  is  proportional  to  the 
square  of  the  combined  amplitude,  we  have 
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/gc-  J a2netdt 

1  -772 
j  772 

=  —  J  [cos(£Z?9  +  kdx  -  co t)  +  cos(kd2  -  co t)]2dt 
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Figure  1.  Schematic  Diagram  of  an  Interferometer 

Using  the  identities  cosv  cosy  =  ^[cos (jc+y)  +  cos(v-y)]  and  cos2jc  =  1/2[1+cos2jc],  we  get 

j  7-/2 

/  oc  —  J  [  l  +  l  /  2  cos(2£2?9  +  2 kdx  -  2cot)  +  1/2  cos(2kd2  -  2cot)  + 

T  _T/2 

cos(kBQ  +  k[dx  +d2\-  2co t)  +  cos(kBQ  +  k[dx  -  d2  ])\dt 

The  second,  third,  and  fourth  terms,  which  are  cosines  whose  angles  contain  an  co?  term, 
will  integrate  to  zero,  leaving 
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-j  T/2 

I  gc  —  J[l  +  cos(kB9  +  k[d1 -d2])]dt 

^  -T/2 


or  simply 


I  oc  1  +  cos(kB9  +  k[d,  -  d2 ]) . 


If  the  flux,  or  energy  per  time  per  area,  at  both  apertures  is  F,  and  if  the  apertures  each 
have  equal  area  A,  then  the  time-averaged  power  from  the  point  source  is 

/  =  2AF[l  +  cos(F89  +  8)]  ,  (1) 

where  8  =  k(d,  -d2)  . 


Equation  1  generates  a  fringe  pattern  that  oscillates  between  9  and  4AF,  where  AF 
is  the  power  collected  by  one  of  the  apertures.  We  could  adjust  the  relative  delay  8  to  a 
value  -FB9  ±  n2r,  where  n  =  9,  1,  2,  ...  ,  to  obtain  the  fringe  at  its  maximum,  represent¬ 
ing  total  constructive  interference.  Total  destructive  interference  is  obtained  when  8  =  - 
kBQ  ±  nr,  where  n  =  1,  2, .... 

Now  suppose  that  a  second  distant  quasi-monochromatic  source  is  illuminated  at 
an  angle  9'  from  the  bore  sight  of  the  interferometer.  Assume  that  the  sources  are 
incoherent,  meaning  that  the  fields  from  the  two  sources  are  uncorrelated,  and 
consequently,  we  do  not  have  to  consider  interference  effects  between  the  sources.  '  Each 
source  will  therefore  contribute  its  own  fringe  pattern  in  the  interferometer,  and  the 
fringes  of  this  second  source  will  be  offset  from  the  fringes  of  the  first  because  9'  ^  9  . 
The  overall  fringe  pattern  will  simply  be  the  sum  of  the  fringes  of  each  individual  source. 
Thus,  /  will  never  reach  a  maximum  of  4  times  the  power  collected  by  one  aperture,  nor 
will  it  fall  to  a  minimum  of  9,  so  the  fringe  pattern  from  the  pair  of  offset  point  sources  is 
not  as  distinct  as  it  would  be  if  the  pair  of  sources  were  collocated  and  together  formed  a 
point  source.  The  “distinctness”  of  the  fringes,  referred  to  as  the  visibility,  is  defined  as 

_  Anax  A n i n 

1  rnax  An  in 


A  point  source  is  characterized  by  V  =  1,  and  a  pair  of  separated  point  sources  (or  any 
extended  object,  for  that  matter),  is  characterized  by  V  <  1. 


1"  A  given  point  source  is  coherent  with  itself — this  is  what  gives  rise  to  the  interferometer’s  interference 
effects  summarized  in  Equation  1.  The  subject  of  coherence  will  be  covered  in  more  depth  in  the  next 
chapter. 
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The  response  of  an  amplitude  interferometer  can  be  used  to  do  more  than  simply 
distinguish  between  point  and  extended  sources — it  can  be  used  to  measure  the  angular 
size  of  objects  and  can  actually  synthesize  the  image  of  an  object.  These  concepts  are 
developed  in  the  following  sections. 


B.  RESPONSE  TO  AN  EXTENDED  SOURCE 

Consider  a  one-dimensional  extended  source  of  incoherent  light.  The  flux  of  the 
extended  source  at  the  apertures  is  F  =  F(0),  and  Equation  1  is  valid  for  any  given  point¬ 
like  element  on  the  source: 

/(0)  J0  =  2  AF(0)[l  +  cos(A/?0  +  8)]  <70  . 

Because  the  source  is  incoherent,  each  F(0)r/0  on  the  source  is  uncorrelated  to  any  other, 
and  the  interferometer’s  response  from  the  entire  source  is  an  integral  over  7(0)r/0: 

7  =  J  7(0)  <70  =  2Aj  F(0)  [l  -  cos(kB0  +  5)]  d6  . 

Using  the  identity  cos(^l  ±B)  =  cos  A  cos  B  +  sin  A  sin  B  ,  we  can  rewrite  the  response  as 
7  =  2 A  F(0)  dQ  +  cos  5  J  F(6)  cos  kBO  dQ  -  sin  5  J  F(0)  sin  kBQ  i/d]  , 


or 

7  =  2^4  [  J  F(0)  dO  +  Re{  em  J  F(Q)e,kmdQ  } 


(2) 


where  Re{  }  stands  for  the  real  part  of  the  expression  in  parentheses.  Equation  2  can  be 
normalized  to  Io,  the  total  power  received  by  the  two  antennas: 


7  /  70  =  1  +  Re< 


,5  J  F(Q)eikmdQ 
J  F(Q)dQ 


(3) 


The  complex  visibility  is  defined  as 

F(0)  eikmdQ 

-r -  ,  (4) 

J  F(Q)de 

which  is  the  normalized  Fourier  transform  of  the  flux  distribution  F(0).  This  definition 
can  also  be  recast  in  terms  of  |T(F)|  and  0(7?),  the  modulus  and  the  phase  of  the 
complex  visibility: 

r(B)  =  |r(F)|  emH)  .  (5) 

The  modulus,  or  amplitude,  which  is  bounded  by  0  and  1,  can  be  written  as 
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r  (5) 


1 


[(] 

1"  5(9)  cos  £59  dQ 

Mj 

1"  5(0)  sin  £59  d@) 

f? 

1 

5(0) 

dQ 

(6) 


The  phase  gives  the  relationship  between  the  real  and  complex  parts  of  the  transform. 

Substituting  Equations  4  and  5  into  Equation  3,  we  get  the  normalized  interfer¬ 
ometer  response  to  be 

I/I0  =  l  +  i?e{e'6  |r(5)|eW(a)}  , 


or  simply 

I/I0  =  l  +  |r(5)|cos(8  +  0>(5))  . 


(7) 


The  interferometer  fringe  response  given  in  Equation  7  is  pictured  in  Figure  2.  The  phase 
0(5)  is  the  offset  of  the  fringe  modulation  from  8  =  0  (which  is  the  interferometer’s 
optical  path  difference  of  0).  The  quantity  |r(5)|  is  the  amplitude  of  the  fringe 
modulation.  In  fact,  the  amplitude  can  be  written  as 


E(5) 


Anax  /„i„i  _  y 

Anax  Anil! 


(8) 


The  modulus  of  the  Fourier  transform  of  the  flux  distribution  is  the  visibility  of  the 
interferometer  response  as  defined  above. 


C.  THE  AMPLITUDE  INTERFEROMETER  AS  AN  IMAGING  DEVICE 

Equation  4  shows  that  a  source’s  flux  distribution — or  “brightness”  distribution  as 
it  is  called  in  astronomy — is  the  inverse  Fourier  transform  of  the  complex  visibility  of  a 
distant,  incoherent  source  that  subtends  small  angles.  This  relationship  is  known  in  optics 
as  the  van  Cittert-Zemike  theorem  [1].  Since  the  complex  visibility  is  measured  by  an 
interferometer,  the  interferometer’s  response  can  be  inverted  to  obtain  a  brightness  image 
of  a  source.  Astronomers  have  coined  this  process  synthesis  imaging  [2,  4]. 
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Figure  2.  Response  of  an  Amplitude  Interferometer  to  an  Extended  Source 

The  difficulty  level  of  synthesis  imaging  depends  on  the  wavelength  of  the  source 
radiation.  This  is  primarily  due  to  measuring  the  phase  ®,  which  is  affected  by 
atmospheric  turbulence  and  instrumental  instabilities.  Although  such  measurements  are 
possible  in  the  radio  frequency  (RF)  spectrum  where  wavelengths  are  macroscopic,  it  is 
quite  a  different  story  at  optical  wavelengths  that  are  a  fraction  of  a  micron.  The  quantity 
V,  on  the  other  hand,  is  more  readily  measured  at  a  variety  of  wavelengths  from  optical 
through  RF  wavelengths.  One  could  imagine  tracking  and  measuring  the  peaks  and 
troughs  in  V  in  Figure  2,  even  as  errors  in  the  optical-path  difference  in  the  interferometer 
make  the  measurement  of  ®  impossible.  At  small  wavelengths,  however,  even 
measurements  of  V  are  challenging. 

So  it  is  not  surprising  that  synthesis  imaging  of  the  sort  described  above  is 
generally  the  domain  of  radio  astronomy,  where  antenna  arrays,  such  as  the  National 
Radio  Astronomy  Observatory  (NRAO)  in  Socorro,  New  Mexico,  have  been  used  to  map 
cosmic  radio  emissions  for  decades. £  Synthesis  imaging  at  optical  wavelengths  has  only 
recently  been  undertaken  [2].  (Note  that  synthesis  imaging  at  optical  wavelengths  should 


I-  Background  information  on  the  NRAO  facility  and  several  extraordinary  radio  images  are  available 
through  the  NRAO  website  (http://www.nrao.edu/). 
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not  be  confused  with  standard  optical  imaging,  where  a  single  aperture  and  lens  or  lens 
group  focuses  light  on  an  image  plane.) 

Although  full  synthesis  imaging  is  exceedingly  difficult  at  optical  wavelengths,  it 
is  possible  to  measure  V  as  discussed  above.  It  has  long  been  known  that  the  angular  size 
of  an  object  can  be  extracted  from  V  in  a  model-dependent  manner  when  the  object’s 
shape  is  known.  This  will  be  shown  in  the  next  section.  It  is  this  technique  that  enabled 
Michelson  and  Pease  and  others  to  measure  the  diameter  of  stars  a  century  ago. 

D.  VISIBILITY  AND  SOURCE  MORPHOLOGY 


1.  Point  Source 

Consider  a  point  source  located  at  0O: 

^point(0)  =  ^5(0-0o)  • 

Replacing  Fpoint  in  Equation  4  gives 

T(B)  =  em%  , 

which,  when  rewritten  in  the  form  of  Equation  5,  gives  |r(5)|  =  1  and  O  =  kBQo-  This 
simply  gives  the  same  result  as  in  Equation  1.  The  visibility  V,  given  by  |E(F?)|,  is  1  for  all 
values  of  B ,  so  it  is  said  that  the  source  is  unresolved  by  the  interferometer  (i.e.,  no 
physical  extent  is  detected). 


2.  Binary  System  of  Point  Sources 

Consider  a  binary  system  of  point  sources  separated  by  angle  a  and  centered  at 
0o: 

W@)  =  A[  8(8  -  (0,  +  a/2))  +  6(8  -  (90  -  a/2))  ]  . 


Replacing  Fbmary  in  Equation  4  gives 

r(s)  =  4k 


ikB(Q0  +  a/ 2)  ^ikB(Q0- a/2) 


=  e 


ikBQn 


ikBa/  -ikBa/ 
,  /2+e  /2 


ikBQ0  „„„ 

=  e  0  cos 


(kBa4 


When  rewritten  in  the  form  of  Equation  5,  we  get  the  visibility  T(F)|  = 
the  phase  <D  =  kBQ0-  The  visibility  of  the  binary  system,  shown  as  a  solic 


cos 


kBa / 


and 


line  in  Figure  3, 


reaches  its  first  minimum  at  a  baseline  of  B  =  a/2  a.  It  is  said  that  this  source  is  resolved 
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by  the  interferometer  because  the  visibility  changes  with  baseline,  unlike  the  single  point 
source  case  where  the  visibility  is  1  for  all  baselines. 


B 


Figure  3.  Visibility  versus  Baseline  for  (1)  a  Binary  System  of  Point  Sources  Separated  by 
Angle  a,  (2)  a  Flat  Source  that  Subtends  an  Angle  a,  and  (3)  a  Circular  Source  (Disk)  that 

Subtends  an  Angle  a 


3.  Flat  Source  Distribution 

Consider  a  flat  source  distribution  with  an  angular  size  a  and  centered  at  Qo'- 

^,(0)  =  —  when  e0--^<e<e0  +  -^ 

a  2  2  • 


=  0 


otherwise 


Replacing  Fflat  in  Equation  4  gives 


T  (5)  =  —!—(< 
ikBa  v 


ikB(Q0  +a/2)  _  ikB(Q0-a/2)\ 

t-  I 


=  e 


ikBQa 


1 


=  e 


ikBa ' 

,km„  sin(CSa/2) 


ikBa/  -ikBa/ 

e  /2  -  e  /2 


(kB  a/ 2) 


When  rewritten  in  the  form  of  Equation  5,  we  get  the  visibility  r(Z?)  = 
and  the  phase  <E  =  kBQo-  The  visibility  of  the  flat  system,  shown  as  a 


sin(kBa/2) 

(kba  /  2) 
dashed  line  in 
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Figure  3,  reaches  its  first  minimum  at  a  baseline  of  B  =  XI  a.  Again,  this  source  is 
resolved  by  the  interferometer  because  the  visibility  changes  with  baseline. 


4.  Circular  Source  Distribution 


Consider  a  uniformly  bright  circular  source  distribution  (a  disk)  with  an  angular 
size  a,  centered  on  the  interferometer’s  bore  sight  (0o  =  0),  as  pictured  in  Figure  4.  The 
denominator  of  the  visibility  integral  in  Equation  4,  which  integrates  Fdisk(0)  over  the 
extent  of  the  source,  is  simply 


a/2 


Jf(0)</0  =  J/?(0>/0 


-a/2 


where  h  is  pictured  in  the  figure.  It  is  possible  to  change  variables  from  0  to  x  as  shown  in 
the  figure  so  that  F(0)  — >  Fdisk(x)  =  acosx  ,  0  =  (a/2)sinx,  and  dQ  =  (a / 2) cos xdx . 
The  limits  on  the  integral  also  change:  0  =  ±  a/2  — >  x  =  ±  ti/2  .  Thus, 


V-2  a2 


Jf(0)g?0  =  J  — cos 2  xdx  = 

-  k/2  ^ 


a  n 


(9) 


The  numerator  of  Equation  4  is  then 

jF(Q)e‘kmdQ  =  J  —cos 2  x  e‘kBa/2  dx  = 


— 7t/2 


xca 

kB 


./,  (7v/?a/  2 ) 


(10) 


Combining  Equations  (9)  and  (10)  gives 

JdiBa/2) 

(kBoJl) 


When  rewritten  in  the  form  of  Equation  5,  we  get  the  visibility  for  a  uniformly  bright  disk 
J , (kBa / 2 )  j,  w|1jc|1  js  s)lown  as  a  dotted  line  in  Figure  3.  The  disk’s 


as 


|F(2?)|  = 


(kba/2) 

visibility  reaches  its  first  minimum  at  a  baseline  of  B  =  1.22X1 a.  Again,  this  source  is 
resolved  by  the  interferometer  because  the  visibility  changes  with  baseline. 
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Figure  4.  Geometry  for  a  Circular  (Disk-Shaped)  Source  that  Subtends  an  Angle  a 
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III.  INTENSITY  INTERFEROMETRY 


As  discussed  above,  measurements  with  an  amplitude  interferometer  can  be 
difficult.  R.  Hanbury-Brown  and  R.Q.  Twiss  made  the  remarkable  realization  that  the 
modulus  of  the  complex  visibility  of  a  source  can  be  measured  by  correlating  the 
intensities  of  two  detectors  separated  in  space  [3].  As  will  be  seen,  there  are  important 
advantages,  as  well  as  some  drawbacks,  to  using  this  approach. 

Figure  5  shows  a  schematic  of  an  intensity  interferometer.  The  outputs  of  two 
square-law  detectors^  are  fed  into  low-pass  filters  and  then  correlated.  The  correlation 
function  computed  at  the  final  stage  is  effectively 

C(x)  =  ±-\la{t)Ib(t  +  x)dt  ,  (11) 

^  0 

where  Ia  and  //,  are  the  filtered  intensities,  and  the  maximum  of  the  correlation  is  at  x  =  0, 
which  corresponds  to  the  optical  path  difference  of  0.  The  correlation  is  simply  the  time 
average  of  la  and  Ih  over  time  T. 

Before  deriving  the  relationship  between  the  correlation  of  intensity  fluctuations 
and  the  visibility  of  a  source,  we  consider  a  simple  model  for  an  incoherent  source,  which 
is  important  to  the  derivation  of  the  Hanbury  Brown  and  Twiss  (HBT)  effect. 

A.  MODEL  FOR  AN  INCOHERENT  (THERMAL)  SOURCE 

The  response  of  an  amplitude  interferometer  is  typically  derived  by  considering  a 
monochromatic,  incoherent  source.  A  monochromatic  source  is  by  definition  coherent, 
however,  so  the  incoherence  in  the  derivation  of  the  amplitude  interferometer  is  included 
in  an  ad  hoc  manner  by  ignoring  the  correlations  between  different  parts  of  the  source.  A 
source  with  a  finite  bandwidth  is  needed  to  understand  the  HBT  effect. 

A  wave  can  be  defined  as  coherent  if  it  is  possible  to  predict  its  phase  over 
arbitrarily  long  times.  If  it  is  only  possible  to  predict  a  wave’s  phase  over  time  ~  xc,  it  is 
said  that  the  wave  is  partially  coherent  with  a  coherence  time  xc. 


§ 


A  square-law  detector  is  one  whose  output  is  proportional  to  the  intensity  of  the  input  field. 
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Figure  5.  Elements  of  an  Intensity  Interferometer 

Finite  coherence  time  implies  finite  bandwidth.  To  see  this,  consider  a  constant- 
amplitude  wave  train,  F(t),  that  has  phase  stability  for  time  periods  exactly  rc.  The 
spectral  density  of  this  wave  train  is  based  on  its  basic  part,  namely  a  harmonic  wave 
with  frequency  coo  and  amplitude /o: 

Fit)  =  foe  iw°‘  \t\<^ 

F{t )  =  0  kl  >  —  • 

M  2 


The  spectral  density, /(co),  of  Fit)  is  given  by  its  Fourier  transform: 

Xc/2 

fit o)  =  J  eimF(t)  dt=  |  /Oe'(co_t°o)/  dt 


-Tc/2 


=  foe 


sin(co-co0)Tc  /  2 
(co-co0)rc/2 


2 

The  intensity  distribution  is  then/(co)L  which  has  its  first  zero  when  the  argument  of  sine 
is  equal  to  7T,  corresponding  to  the  width 

(co-co0)tc  /2  =  n 
or  Av  =  1  /  xc  . 


Thus,  the  bandwidth  is  related  to  the  finite  coherence  time. 

For  the  purposes  of  this  study,  an  incoherent  wave  train  can  be  constructed  by 
using  a  monochromatic  wave  with  random  phase  jumps,  uniformly  distributed  between  0 
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and  27T,  that  occur  on  average  every  xcM  The  actual  distribution  in  time  for  the  phase 
jumps  is  not  important;  a  flat  or  Gaussian  distribution  with  a  mean  of  xc,  or  even  a  phase 
change  every  xc,  would  work  to  develop  the  HBT  concept. 

B.  DERIVATION  OF  THE  HBT  EFFECT 

Figure  6  is  a  schematic  diagram  of  two  identical  detectors,  a  and  b,  viewing  a 
source  consisting  of  many  points  such  as  m  and  n.  The  detectors,  which  measure 
intensity,  are  read  out  as  pictured  in  Figure  5.  The  field  due  to  point  m  at  detector  a  is 
given  by 

K,  =  A„,  cos(A-/r  -  ov  +  <k,  )  , 

where  the  phase  4>m  is  a  time-dependent  quantity  ( <j)m  =  <j),„(0)  that  changes  with  a  period 

that  is,  on  average,  equal  to  the  source’s  coherence  time  xc.  The  phase  changes  are 
randomly  distributed  between  0  and  2ir. 


Source 


Detector  a  Detector  b 

Figure  6.  Schematic  Diagram  of  an  Intensity  Interferometer  Viewing  a  Source 

The  field  from  the  whole  source  at  a  is  then 


V  This  construct,  though  an  oversimplification,  does  access  an  important  characteristic  of  a  thermal 
source,  namely  that  the  phase  of  an  emitting  atom  is  interrupted  because  of  atomic  collisions.  Although 
there  are  other  effects  in  thermal  emissions,  such  as  the  damping  of  wave  trains  as  atoms  lose  energy 
during  emission  and  the  Doppler-broadening  of  lines  [1],  these  are  not  necessary  to  understanding 
HBT. 


HI-3 


A“  =IX  cos(W2“ -orf  +  <|>J 

m 


and  the  intensity,  Ia,  is 

l .=(^“)2=  Z4,  «»(**; -«*  +  ♦„) 


Z4,  cost /  /?;;  01/ 1  o.: ) 


4zi4,4,W*(*;-s;)+i  -*.)+cos(*(tf; +o-2«x+*. +0] 

^  m  n 

where  we  have  used  the  identity  cosv  cosy  =  V2[cos(x+y)  +  cos(v-y)].  Assuming  the 
detector’s  output  is  low-pass  filtered,  we  are  left  with 


I  =-YYA  A  cosfrcT -Ra)  +  6  -d>  ) 

a  r\  m  n  \  V  m  ns  >  m  t  n ) 

^  m  n 

=  ~Y  A2  +  ~Y  A  A  cos[k(Ra  -Ra)  +  6  -(b) 

^  /  -<  /w  ^  m  n  \  \  m  ns  t  m  \  n ) 


(12) 


We  are  interested  in  computing  Equation  11,  which  is  the  time  average  of  the  intensities 
of  detectors  a  and  b.  The  time  average  is  much  longer  than  the  coherence  time  of  the 
source.  If  we  denote  the  time  average  of  the  intensities  at  r  =  0  with  angle  brackets,  we 
have 


<IaIb>=\{  terms  linear  m  cosine 

^  \  m  m' 

+  Z  Z  4, ,4,4/4/ coslHK-K )  +  *„- 


m^n  m'^n' 

•cos[*(J&-#)  +  4>b.-(M 


(13) 


The  terms  linear  in  cosine  will  time  average  to  0  because  the  time  average  is  much  longer 
than  the  coherence  time,  hence  jumps  in  phase  (<j))  will  cause  the  cosine  to  average  to  0. 
The  last  term  in  Equation  13  can  be  rewritten  as: 


X  X  AmAnAm,An,  cos[£(i?“  - iC)  +  (L  - <U  cos[*(J&-#)  +  *B.-(U 


\m^n  m  ^ n 


Z  Z  4, 4,4/4/  {kos[*-(A,:  -  r:  -  R‘.  +  Rbn.)  +  +.  -  <|>„  - +  4v] 


.  (14) 


\m^n  m  ^ n 


+ -cosc^i?:  -  r: + ri  -  k.) + <►„  -  i + *m.  -  i.] 


The  first  term  on  the  right-hand  side  of  Equation  14  survives  the  time  average  only  if  m  = 
tri  and  n  =  ri,  or  else  the  phases  (c|)  terms)  remain,  which,  over  time,  cause  the  cosine  to 
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average  to  0.  Similarly,  the  second  term  survives  only  if  m  =  n'  and  n  =  in'.  Thus 
Equation  14  reduces  to 


Z  Zaaaa  cos[k(x: -r:)+ -  cUcos^r*,,  -  K) + $m.  -  cu 

=  Y  AlAl  cos [k(Ram  -  Rbm  -  R°  +  Rbn)] 

m  n  l  v  m  m  n  ns. 1 


Substituting  this  result  into  Equation  13  gives 

(  4 

I  1^  o 

<44  >= 


+7Z-4»4?coS[*-«-^,-<+^)]  .  (15) 


m  n 


\  m  J 

The  time-averaged  power  from  one  detector,  as  given  by  Equation  12,  is  simply 


<4  >=\ZAl’ 


(16) 


because  the  leftover  phases  cause  the  second  term  of  Equation  12  to  time  average  to  0. 
Both  detectors,  a  and  b,  will  measure  the  same  time-averaged  power  as  given  in  Equation 
16,  which  can  be  called  simply  </  >.  We  can  then  normalize  Equation  15: 

.  £ aIa; cos [*(*; - r‘ - k  +  sp] 

a  b  _  |  j  m^n _ 


<I>Z 


1 4 

V  m  J 


Inspection  of  Figure  7  shows  that  Ram  -  Rbm  =  BQm  ,  thus 


<  44  ^  _  Y  +  m*n 

<  I  4 


(  V 

IX 

V  m 


Rewriting  the  limits  on  the  summation  gives 

2X  X44!cos[M(e,-ej] 


<  i,h>  , 

</>2 


(  V2 

ZAi 

V  m  J 


+  - 


(  V2 

IX 

V  m  J 


(17) 
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Source  point  m 


Figure  7.  Path  Length  Difference  for  Radiation  Arriving  at  Detectors  a  and  b.  B  is  the 
baseline  of  the  two  detectors  and  6  is  the  angle  of  source  point  mfrom  the 

interferometer’s  bore  sight. 


Consider  the  second  term  on  the  right-hand  side  of  Equation  17.  Assume  that  the  source 
to  be  imaged  consists  of  N  equal-strength  point  sources  and  that  the  interferometer  is 
capable  of  resolving  the  points  within  the  source  (i.e.,  )JB  is  small  enough  to  do  this). 
The  second  term  in  Equation  17  then  reduces  to 


1 4 


m 


i 

7  \2  N 

2 


V  m  J 


which  vanishes  in  the  limit  of  a  source  built  of  many  points.  Another  way  to  view  this  is 
that  if  the  interferometer  resolution  is  small  compared  to  the  angular  size  of  the  source, 
and  if  the  changes  in  brightness  are  not  too  dramatic  across  the  source,  then  this  term  can 
be  ignored.  This  term  cannot  be  ignored  if  the  source  to  be  imaged  is  made  of  a  small 
number  of  points,  or  if  it  is  made  of  many  points  but  with  a  few  that  are  much  brighter 
than  the  rest.  Nevertheless,  for  many  interesting  cases — where  complex  sources  are 
viewed  by  interferometers  with  adequate  resolution — this  term  can  be  ignored. 
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The  third  term  in  Equation  17  can  be  rewritten  using  the  definition  for  </>  and  the 
identity  cos  (A  — B)  =  cos  A  cos  B  +  sin  A  sin  B  : 


X4;4>s[«(8.-e,,)] 

m,n 

f  V 

\  m  ) 


1 

<  /  >2 


Z44(coskB0m  cosA7?0„  +  sin kBQm  sin kBQn) 

m,n 


1 

</>2 


Z^* cosAfi0™ 


\  m 


J 


(  3 

2 

+ 

Z  A;„  sin  A7?0m 

V  m  J 

2 

Noting  that  Am  is  just  the  intensity  (or  brightness)  Fm  of  the  point,  and  converting  sums 
into  integrals  by  replacing Fm  — >  F(Q)dQ ,  the  intensity  correlation  in  Equation  17  can  be 
rewritten  as 


<  I  ah  > 
</>2 


=  1  + 


<I>Z 


(j/r(0)cos£60d0)2  +  (jFCOjsin^Od©)2 


A  comparison  of  this  result  to  Equation  6  shows  that  the  intensity  correlation  is  given  by 

<44  >  i  +  ir(i?)|2  (18) 

<i>2  1  1 

The  intensity  correlation  measures  the  modulus  (but  not  the  phase)  of  the  Fourier 
transform  of  the  source  distribution. 


C.  ADVANTAGES  AND  LIMITATIONS  OF  THE  HBT  EFFECT 

HBT  cannot  be  used  to  fully  image  a  complex  source  because  it  measures  only  the 
modulus  (and  not  the  phase)  of  the  complex  visibility.  Knowledge  of  the  phase  is 
necessary  to  invert  the  complex  visibility  and  obtain  the  source’s  brightness  distribution. 
HBT  can  be  used  to  measure,  in  a  model-dependent  fashion,  the  size  of  objects  whose 
shape  is  known.  An  obvious  application  of  this  is  the  measurement  of  the  angular  size  of 
stars.  It  can  also  be  used  to  obtain  the  approximate  size  of  an  object  whose  shape  is  not 
known.  This  can  be  seen  by  examining  Figure  3,  which  shows  the  modulus  of  the 
complex  visibility  for  different  source  shapes.  Although  the  specific  source  shape  will 
determine  the  precise  baseline  at  which  the  magnitude  of  the  visibility  vanishes,  HBT 
measurements  can  give  an  approximate  measurement  of  the  size  of  the  overall  source — 
the  first  minimum  in  visibility  for  a  binary,  flat,  or  circular  source  is  within  a  factor  of  2.5 
of  one  another  as  shown  in  the  figure. 

Within  the  context  of  measuring  the  modulus  of  the  complex  visibility,  there  are 
advantages  and  disadvantages  of  using  HBT  instead  of  amplitude  interferometry.  One 
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advantage  is  that  HBT  is  less  sensitive  to  atmospheric  turbulence  and  device  instabilities. 
In  amplitude  interferometry,  the  phase  of  the  incident  radiation  after  having  taken 
separate  paths  through  the  medium  is  measured  at  two  separate  detectors.  A  fluctuating 
medium  or  device  instabilities  that  induce  path  length  errors  of  order  a  wavelength  would 
seriously  affect  the  accuracy  of  the  measurement.  The  HBT  interferometer  need  not 
preserve  the  stability  at  the  level  of  the  incident  radiation  wavelength.  Instead,  stability 
only  on  the  level  of  the  coherence  time  is  needed,  hence  the  interferometer  is  less 
sensitive  to  atmospheric  effects  and  device  instabilities.  HBT’s  disadvantage  is  that  it  has 
diminished  sensitivity  for  dim  sources  because  the  signal-to-noise  ratio  at  the  HBT 
correlator  goes  as  (PS/PN)2,  where  Ps  and  PN  are  the  signal  and  noise  power  at  detectors. 
In  an  amplitude  interferometer,  the  overall  signal-to-noise  ratio  goes  as  Ps/Pn- 

D.  EXTENSION  OF  HBT 

HBT  measurements  cannot  image  a  complicated  source  because  they  do  not 
provide  the  phase  of  the  Fourier  transform  of  a  source’s  brightness  distribution;  however, 
methods  have  been  developed  to  extract  the  phase  information  in  intensity  correlations. 
Gamo  proposed  the  superposition  of  a  coherent  beam  on  the  measurement  of  an  unknown 
source  [5].#  In  another  approach,  Mehta  proposed  that  a  source  with  known  but  finite 
coherence  be  used  to  extract  the  phase  data  [6].  This  approach  was  demonstrated  in  the 
lab  [7].  Although  resources  preclude  developing  these  ideas  in  detail  here,  we  offer  a 
brief  summary  that  follows  [6].  The  visibility  of  an  unknown  source  is  measured  to  be  V 
using  standard  HBT  measurements.  A  reference  source  with  known  visibility  and  phase, 
Vo  and  ®o,  is  superimposed  on  the  unknown  source,  and  the  visibility  of  the  combined 
source,  V ,  is  also  measured  using  standard  HBT  techniques.  The  visibility  of  the 
combined  source  is  given  by  V 2  =  V 2  +V02  +2 W0  cos(®  -®0),  where  ®  is  the  phase  of 
the  unknown  source.  Since  all  of  the  quantities  except  ®  in  the  expression  for  V2  are 
either  known  or  measured,  ®  -  ®o  can  be  solved  to  within  a  factor  of  -1.  The  sign 
ambiguity  results  in  an  inverted  and  uninverted  image  [7].  This  process  is  analogous  to 
holography.** 


*  Gamo  further  proposed  extracting  the  lost  phase  information  by  using  triple  correlations,  or 
correlations  of  three  detector  elements.  See  [5]  for  details. 

T.  D.  Beard,  private  communication. 
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IV.  SUMMARY 


An  amplitude  interferometer  measures  the  modulus  and  phase  of  the  Fourier 
transform  of  the  brightness  distribution  of  an  incoherent  source,  thus  the  response  of  the 
amplitude  interferometer  can  be  inverted  to  obtain  an  image  of  the  source.  This  imaging 
process  is  dependent  on  the  ability  of  the  device  to  accurately  track  the  phase  of  the 
incident  radiation  field.  A  traditional  intensity  interferometer  measures  only  the  modulus 
of  the  Fourier  transform  of  a  source’s  brightness  distribution.  It  cannot  fully  reconstruct  a 
complex  source  because  the  device  does  not  measure  phase  information  of  the  Fourier 
transform.  This  type  of  device  is  nevertheless  able  to  approximate  a  source’s  size,  and,  in 
the  case  where  the  source’s  shape  is  known,  a  precise  size  estimate  can  be  extracted  in  a 
model-dependent  manner. 

Extensions  to  the  intensity  interferometer  enable  this  type  of  device  to  fully  image 
a  source.  These  approaches  involve  the  superposition  of  a  known  source  onto  the  device 
or  processing  triple  correlations  among  interferometer  elements  (as  opposed  to  double 
correlations). 

It  is  appropriate  to  use  an  amplitude  interferometer  when  the  phase  of  radiation 
can  be  measured.  A  source  can  be  fully  imaged  with  this  type  of  interferometer,  and  the 
signal-to  noise  ratio  for  dim  sources  is  higher  in  this  device  than  in  intensity 
interferometers.  If  tracking  the  phase  of  radiation  is  difficult,  as  when  dealing  with  small 
wavelengths  or  atmospheric  or  device  instabilities  that  cause  phase  measurement  errors, 
then  intensity  interferometers  are  appropriate.  Although  the  signal-to-noise  ratio  of  an 
intensity  interferometer  suffers  when  the  source  is  dim,  it  is  more  robust  in  a  turbulent 
medium.  Traditional  intensity  interferometers  are  limited  as  imaging  devices;  however, 
extensions  to  the  intensity  interferometer  allow  for  sources  to  be  fully  imaged. 
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